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SUMMARY 


The development of the equations of dynamic equilibrium for a lifting sur- 
face from Lagrange's equation is reviewed and restated for general exponential 
growing and decaying oscillatory motion. Aerodynamic forces for this motion are 
obtained from the three-dimensional supersonic kernel function that is newly gen- 
eralized to complex reduced frequencies. An existing computer program for the 
kernel-function method was modified to solve the equations. Illustrative calcu- 
lations were made for two flutter models at supersonic Mach numbers. Preflutter 
and postflutter motion isodecrement curves were obtained. Variations of motion 
decrement with Mach number and with air density were much milder than those pre- 
dicted from conventional curves of velocity plotted against damping required 
(the V-g type) using simple harmonic aerodynamics. This type of analysis can 
be used to predict preflutter behavior during flutter testing and to predict 
postflutter behavior for use in the design of flutter suppression systems. 


INTRODUCTION 

The phenomenon analyzed here is the preflutter (decaying or subcritical) 
and postflutter (growing) oscillations of an aerodynamic surface in an airstream, 
which is a nonconservative system. For the usually analyzed condition of neu- 
trally stable flutter motion, the equations of motion express a condition of 
zero net exchange of energy between the fluttering aerodynamic surface and the 
passing air stream. In contrast, for growing and decaying motion, a net contin- 
uing exchange of energy occurs. The approach using Lagrange's equation is 
reviewed and applied for this motion so as to show the way the associated aero- 
dynamic and structural damping forces are accommodated. This approach is done 
in accord with the original application in reference 1 and with the reassurance 
given by reference 2. 

Hassig in reference 3 described what is called a true-damping solution for 
growing or decaying motion and presented sample results obtained by using aero- 
dynamic forces for simple harmonic motion. Those results were obtained by three 
solutions called the p-method, the p-k method, and the k-method. Each of 
these solution methods gives near flutter results that are approximate because 
of the simple harmonic aerodynamics employed. Of the three, the p-k method is 
believed to give the nearest approximation to what would be obtained with the 
motion-matched aerodynamic forces, as developed and used in the present paper. 

A few analyses of near-flutter motion with motion-matched aerodynamic 
forces have been reported. Dugundji, Dowell, and Perkin include in references 4 
and 5 analyses and results for the unstable motion of panels at speeds greater 
than the flutter speed by permitting a complex wave speed for the traveling 
waves of panel deflection. Morino presents in reference 6 the development of 
an aerodynamic potential by using a Green's function for timewise arbitrary 
motions. Some applications of this method given in reference 7 show generalized 



aerodynamic coefficients for growing, steady, and decaying oscillations of a 
rectangular wing at supersonic speeds and some chordwise distributions of pres- 
sure coefficient for growing oscillations of a wing -body- tail configuration 
pitching in a subsonic flow. As part of its overall development based on the 
Laplace transform technique, reference 8 gives certain functions for general expo- 
nential motions in two-dimensional flow. These are the generalized Theodorsen 
circulation function for two-dimensional incompressible flow and the velocity 
potential and associated lift and moment functions for two-dimensional super- 
sonic flew. 

In the present paper the aerodynamic forces for the growing and decaying 
oscillatory motion are obtained by generalizing to complex frequency the super- 
sonic three-dimensional velocity potential of reference 9 and the kernel func- 
tion of reference 10. The analysis is illustrated by applying at supersonic 
speeds to two flutter-tested configurations the kernel-function aerodynamic 
forces in a computer program modified from that of reference 11. This type of 
analysis has applications to the prediction of motion decrements during wind- 
tunnel and flight flutter testing and to the design of automatic controls for 
flutter suppression. 

The appendix reviews the relation between the structural damping coeffi- 
cient and the logarithmic decrement of motion in the presence of solid-friction 
structural damping (mechanical hysteresis) as used in the analysis. 


SYMBOLS 

a real part of p (see eq. (2)) 

bjj half of mean geometric chord 

bg half of root chord (see fig. 1) 

c speed of sound 

g structural damping coefficient 

h shape of modal deflection (see eq. (4)) 

i unit of imaginaries, 

K supersonic kernel function 

k reduced frequency, wbg/V 

1 reference length 

M Mach number 

m mass per unit area 

m^ generalized mass of mode i (see eqs. (20)) 
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exponential coefficient of motion in ePt 
lifting pressure, positive with z 
generalized force (see eq. (5)) 

element of generalized aerodynamic-force matrix (see eq. (14)) 
generalized coordinate of motion 
amplitude of q (see eq. (4)) 
area of lifting surface 

ratio of semispan to root chord (see fig. 1 ) 

kinetic energy of structure 

time 

potential energy of structure 
unit function (see eq. (29)) 
free-stream speed 
downwash on lifting surface 

time and space varying functions, respectively, in w (see eq. (24)) 

right-hand Cartesian coordinates, x positive to rear, z positive 
upward 

= X - ? 

= y - n 

z-coordinate of deflection shape 
= \/m2 - 1 

logarithmic decrement of motion (see eqs. (3), (A3), and (A7)) 
virtual displacement and virtual work, respectively 
limits of Ti-integration (see eq. (25)) 
a 

= 1 - i — (see eq. (2)) 

0) 

mass ratio 
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t 
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dummy variable for x and y, respectively, representing locations 
of sources and doublets 

air density 

effective spring stiffness of structure (see eq. (7)) 
variable of integration in K 
limits of integration in K 

times of origination of disturbances that affect point (x,y) at time 
t (see eq. (25)) 

perturbation velocity potential on upper surface 
potential function for doublet located at origin 

= — 1 - i - (see eq. (22)) 


^0) 


■n 


w circular frequency of motion 

w = o)/v32 

Wjj base or reference frequency 

Wi,a3j natural frequencies of modes i and j, respectively 
Subscripts; 

aero contribution fron aerodynamic forces 

b base or reference mode 

i,j modes 

sf solid friction (mechanical hysteresis) 

Mathematical notation; 

(*) first time derivative, 3( )/3t 

(■■) second time derivative, 9^( )/3t^ 
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ANALYSIS 


General Exponential Motion 

For the purpose of this paper it is sufficient to analyze the perturbations 
of an initially plane surface. Figure 1 shows the lifting surface and the coor- 
dinate system. Essentially all findings apply to nonplanar and multiple lifting 
surfaces as well, but the procedure is illustrated without the complications of 
nonplanar and multiple surfaces. The surface analyzed has its timewise mean 
location in the plane z = 0. The z-coordinate of the instantaneous perturba- 
tion away from the mean location of the surface is 

z(x,y,t) = Z(x,y)ePt (1) 

where the space variation Z(x,y) can be complex. The exponential coefficient 
p is of the form 


p = a + ico = ito - i — j = io)9 (2) 

where i = is the unit of imaginaries, w is the circular frequency, and 

a is the exponential time rate of growth of motion. For a/u « 1 


fl + i — + 0(62)] 

/ . M 

" i(i) 1 + 1 — 

(3) 

V 2tt / 

\ 2tt/ 



where 6 is the logarithm of the ratio of the amplitudes at the beginning and 
end of one complete cycle and is positive for a decaying oscillation. It is 
commonly called the log decrement. (See appendix.) 


Equations of Equilibrium 

For the general unsteady aerodynamic forces considered here, direct solu- 
tion of the equations of motion for flutter is usually not feasible. Accord- 
ingly the flutter deflection is approximated in the usual way by a finite modal 
series 


z(x,y,t) -x ^ qi(t) hi(x,y) 
i 

qi(t) = qiePt 


(4) 
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where the modal deflections hi(x,y) are real functions as specified in refer- 
ence 2, and is the complex amplitude of qj^(t). 

The equilibrium equations are obtained from Lagrange's equation, which is 
a specialized form of the principle of virtual work 


d - 

fd^\ 

3t 

3U 

dt' 

laqj 

3qi 

3qi 


(i = 1, 2, . . .) (5) 


where T is the kinetic energy, U is the potential energy (elastic only in 
the present use) , and is the generalized force associated with the gener- 

alized displacement q^ and is obtained from the virtual work due to external 
forces (aerodynamic forces and structural damping) . 

The kinetic energy of the structure is 


T = J - m(x,y>^^ <3i hi(x,y)j^ dS (6) 

where S is the area of the lifting surface, and m(x,y) is its mass per unit 
area. The elastic potential energy is represented by 

(x,y)j^ dS (7) 

where o(x,y) represents the effective spring stiffness rate of the elastic 
restoring force. 

Application of Lagrange's equation in the absence of external forces gives 



^ qj I m(x,y) hi(x,y) hj (x,y) dS + ^ qj | a(x,y) hi (x,y) hj (x,y) dS = 
j •'S j •'S 


(8) 


Using the relationships (see eq. (4) ) 
qj = qj (t) = qjePt 

qj = qj (t) = pqjePt = pqj 

qj = qj (t) = P^qj ePt = p2qj = -(toe)2qj 




(9) 
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in equation (8) gives 


I dS — (w j0 ) ^q-j I nihj^ 

j •'S j *4 


hj dS 


( 10 ) 


which is the usual Rayleigh type approximation for replacing the elastic restor- 
ing forces involving a in terms of the mass-inertia terms involving nti)j202. 
Note that in the absence of external forces, a = 0, p = io), and w becomes 
the modal eigenvalue O) j . 

The virtual work is that done by the external forces moving through the 
virtual displacement hi6q;j^; thus, 

6w = 5Waero,i "*■ 

The aerodynamic contribution is expressed initially by 


6w 


aero 


,i ' J ^P(x,y,t) hi(x,y) 6qi(tl 


dS 


( 11 ) 


The lifting pressure distribution Ap is approximated by a modal series paral- 
lel to that for deflection z in equation (4) j thus. 


Ap(x,y,t) ~ ‘^'Pj(x,y) 


( 12 ) 


where I is a user-selected reference length such as 2bQ, and Apj is the 
complex amplitude of the lifting pressure distribution caused by ana associated 
with the modal deflection hj (x,y) . with substitutions of the series of equa- 
tions (4) and (12) 


6w 


aero 


,i = J Y ^Pjj (hi<Sqi) dS 


= Qaero,i'5<3i 
from which is obtained 


Q. 


aero 


,i = f Sr 

•'S j ‘ 


dS 




(13) 


= 2 <3jQ 
j 


ID 


(14) 
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where the generalized aerodynamic-force element 


•ij = J y Apjhi 


This development is valid for real mode shapes h^ as determined in reference 2 
and as discussed in C,17 of reference 12. 

The solid-friction structural damping force is of the same form as in the 
damping term of equation (A5) of the appendix. Thus, 


>5Wsf,i = I a(x,y)- ^ hjqj(-gsf,j)hi6qi 


Use of equation (4) gives 


<5Wsf,i = j ^-ia ^ qjhjgsf,jjhi6qd 


“ Qsf,iA<3i 


where the structural damping coefficient 9sf,j differ among the modes j, 

The contribution to Qi is ' 


2sf,i = -i 23 ^sf,jqj I 


hj dS 


Using the equality developed in equation (10) 


Qsf,i ” 


23 9sf , j I mh^hj 

j •'S 


The equality given by equation (10) applied for the condition of no external 
forces, and the modes and frequencies are thus for the natural undamped or nor- 
mal modes for which h^ and hj are orthogonal with respect to the mass dis- 
tribution; that is. 
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( 20 ) 





(j = i) 


(j ^ i) 


Substitution of these developments into Lagrange's equation (5) gives 


-to202mj^q^ + + igi)miqi - ^ *3jQij = 0 (i = 1, 2, . . .) 

j 


( 21 ) 


In order to put the equilibrium equations into a form for routine solution by a 
complex eigenvalue subroutine divide by 


-0)202(1 + igi) (0)i/0)b) 2nii 


Then the result is (for 6 /it « 1) 


/‘^b\^ 

“^6 + — C + igi)"^ 

V'^i/ 


«3i 


( 0 )b/ 0 )i )2 qj r 

1 — 

mi(l + igi)w2n + i -j ^ ® 


dS = 0 


(i — If 2^ • • •) 


(22) 


where 


^hY 6 \ 

^ ~J 

\o) / V it/ 


A matched solution is obtained when the output value of the log decrement 
from an eigenvalue matches the preselected input value in Apj and in the 

denominator term 1 + i6/TT. (The lifting pressure Apj is a function of 6 
as well as of x, y, M, and k.) For the case of simple harmonic aerodynamics 
(input 6=0), equation (22) reverts to the conventional flutter equation, and 
-6 /it in the eigenvalue is equivalent to and interpretable as the conventionally 
used structural damping coefficient g. 
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Aerodynamic Forces for General Exponential Motion 

The lifting pressure distribution for general exponential motion required 
in equations (15) and (22) can be obtained for supersonic speeds by generaliza- 
tion to ccmplex frequencies of the results of references 9 and 10. 

Velocity potential in supersonic flow.- The velocity potential analysis of 
reference 9 is applied here to the general case of exponentially growing and 
decaying oscillatory motion; that is, with the time variation of downwash 

w(t) = ei“teat = (23) 

where a and u are real, and 


e = 1 


i 


a 

0) 


Where a lifting surface is analyzed as a distribution of source-sink pairs over 
the projection of the lifting surface onto the plane z = 0 (see fig. 1), the 
local strength of the source-sink pairs is proportional to the local downwash 

w(x,y,z=0,t) = W(x,y) w(t) (24) 


(The notation here is mostly parallel to that of ref. 9.) Parallel to equa- 
tion (15) of reference 9, the velocity potential for the upper surface is (in 
the plane z = 0) 


<l’u (x,y,t) 


W(C,n) w(t - Ti) + w(t - T2) 

^ d(2bon) d(2boC) 

-2tt 3 \^(n - m) (ri2 - n) 2bo 


where 


(25) 



ri 2 = y + 


X - C 


M(x - ?)2bo \J(T] - TV|) (ri2 - n) 

T-, = 2bo 


c3^ 


c3 


M(x - ^)2bo \/(ri - ni) (ri2 - n) 

c 32 ^ c3 


2bo 
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I 


6 = \jvfi - 1 

c is the speed of sound, and ^ and ri are dummy variables for x and y, 
respectively. The velocity potential for the lower surface is the negative of 
that for the upper surface. In the integrand the time variation is 

w(t - Ti) t w(t - T 2 ) = + ,a(t-T 2 )^io,(t-T 2 ) 


^ g(a+i03)tg-(a+ia)) (T 2 +T-, )/2j^g(a+io)) (T 2 -Ti )/2 


^ ^-(a+iw) (T 2 -Ti)/ 2 j 


. gla)ete"i“®"’‘o2t>o/o62 2 \l ° ^0 0^ I2JJ 


This time variation is the same as that just before equation (19) of reference 9, 
except that co is replaced by u)0. Thus, for the upper surface 


4>u(x,y,t) 


_gito9t 


IT 



w(5,n) u(xo - e|yo 


)e 


-i6O0M2boxo/c3 


2 


where 


X cos 


|^U)e2bo\/xo2 - a^yo^y {2bon) d(2boC) 
C02 j2bQ\jxQ^ - B^yo^ 


(27) 


XQ = X - 5 

Yo = y - 0 

Acceleration potential and kernel function in supersonic flow .- As shown 
in the preceding section for an exponentially growing (or decaying) oscillatory 


n 



motion with time variation thg velocity potential is the same as that 

for a simple harmonic motion e^t, with o) replaced by wd. 

Given the relationship between velocity potential, acceleration potential, 
and the kernel function shown in equations (8) to (12) of reference 10, the 
same replacement of to 6 for oj,k0 for k, and 0)6 for o) applies throughout. 
This is true for supersonic flow because the integrations back over events orig- 
inating in past time, extend over only a finite time, and are therefore distinct 
and finite. 


Thus the potential for the doublet located at the origin is 
g eio)0te iM%j62bQX ^WLo62bo^x2 - 3^(y^ + 

3 (2bQz) 2bo\^x^ - 3^(y^ + z^) 


(28) 


where to = Equation (28) is a generalization of equation (4) of refer- 

ence 10. The planar (z = 0) supersonic kernel can be used in the form 


K(xo»yo) = 


-2 


V(2bo)2yo2 


g-ik6x - 6|yo|) 


-ik6xo/62 

/Mk0 

cos I yxQ^ - 32yo2, 


\Zxq2 - \32 


ir. 

-ik6|yo|T 


. 

+ 

to 

r 

1 


dT 


where the unit function 


U(XQ - 3|yo|) = 


and 


Tj = — ^ (xQ - M^xq 2 - 32yo2) 

32|yoj ^ ’ 


(xQ - 3 yo > 0) 
(xQ - 3 yo = 0) 


(29) 


Tu = — (xQ + - B^YO^) 

e=|yol ' 

Equation (29) is a generalization of equation (16c) of reference 10. A similar 
generalization applies to the nonplanar kernel for supersonic flow such as that 
presented in reference 13 (with sign corrections as explained in section entitled 
"Results and Discussion") . 
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Method of Solution of Kquations 

Numerical solutions of the equilibrium equations are made for a model with 
known parameters by first selecting values of M, k, and 5 and then solving 
the set of equations (29) for its eigenvalues for each of a selected range 

of values of air density p, usually beginning with P = 0. From each eigen- 
value a frequency o) and output value of decrement 6 are extracted, and asso- 
ciated parameters, such as flutter-speed index or stiffness-altitude parameter, 
can be calculated as desired. Plots of desired results are made for the range 
of p used, with attention given to properly connecting the modal roots. Par- 
ticular attention is given to finding the values of P for which the output 5 
equals the input 6. These values indicate matched conditions. (The important 
matched points are usually those near the flutter crossing with the lowest 
dynamic pressure, called the critical crossing.) For a particular k-value the 
matched-point results may not fall in the desired range of density or mass ratio 
In this case a new k-value can usually be estimated on the basis of a constant 

dynamic pressure and constant frequency.^ This estimate gives k proportional 

to \Jp , which is also proportional to \/T7U. 


RESULTS AND DISCUSSION 

Analytical results have been obtained by using a modified version of the 
computer program in reference 11. This program employs supersonic kernel- 
function aerodynamics. For calculating supersonic kernel-function values for 
complex reduced frequencies, it was convenient to use subroutines based on the 
appendix of reference 13. (In ref. 13 the series summations with the coeffi- 
cients ajn and b,n have incorrect signs that were reversed in the present 
use.) Downwash collocation values and modal deflections for calculating gener- 
alized aerodynamic forces were obtained by interpolating the experimental modal 
deflection data by using a surface-spline computer program based on the analysis 
in reference 14. Certain experimental modal deflection data had sufficient 
irregularity to require smoothing, especially in obtaining streamwise slopes for 
downwash collocation. Accordingly the spline program was modified so that, in 
effect, the modal data were attached to the spline by a uniform set of rather 
stiff springs. 


Clipped-Tip Delta-Wing Model 

Experimental mode-shape, natural-frequency, and mass data for a clipped- 
tip delta-wing model are given in tables I to III of reference 15. The first 
six natural modes were used in the present calculations. (Spot checks for 
modal convergence with nine modes indicate that six modes are sufficient for 
convergence. ) 

Simple harmonic aerodynamics .- For M = 1.3, k = 0.47, g^ = 0, and 6=0 
(simple harmonic aerodynamics), figures 2(a) and 2 (bl display two associated 
results, -5/t\ plotted against flutter-speed index V/bmU2'^ and against P, 
respectively. Figure 2(a) includes a sketch of the wing-model planform. Only 
the mode 2 traces are plotted because this mode gives the critical flutter cross 
ing (lowest V/bj„aJ2\/vi) • In the flutter-speed index ^m and y are as in ref- 
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erence 15; bm is one-half the mean geometric chord, and U is the ratio of 
the wing mass to the mass of air contained in the frustum of a cone that just 
encloses the wing. 

The crossing points at 6 /tt = 0 are the ones that match the input 6=0 
in the aerodynamic forces for simple harmonic motion and give the values of 
V/bjn^ \/y and P at the flutter boundary for this M, k, and other parameters. 
At tile end of the section entitled "Analysis" it was pointed out that, for input 
6=0, -6/ IT is interpretable as the conventional g. Thus the curve in fig- 

ure 2(a) is a V-g type. In some past analyses V-g curves have been inter- 
preted as predicting rates of growth and decay of the oscillatory motion. (See, 
e.g., ref. 16 and the description of the k-method in ref. 3.) The two labeled 
points at 0.05 and -0.05 could accordingly be interpreted to indicate near- 
flutter motion with log decrements of motion of -0.051T and 0.051T, respectively, 
but such an interpretation does not match the input 6 in the aerodynamic 
forces. In figure 2(b) the three points discussed (6 = 0, -6 /tt = 0.05 and 

-0.05) are projected to a plot of P against V/b^0)2 \/y. (The latter scale is 
at the upper right.) These three points are used in subsequent figures. 

Matched aerodynamics for growing and decay ing motion, 6 f 0.- Now con- 
sider the result of using input values -6/n = +0.05. Figures 3(a) and 3(b) 
show the same type of results as those in figures 2(a) and 2(b). The mode 2 
traces for 6=0 are replotted here. Segments of the traces for input 
_6/ir = +0.05 are plotted, and the matched crossing points are indicated. At 
the lower motion-matched point -6/fr = -0.05, the decay rate for the predicted 
motion matches that of the aerodynamics. Similarly at the upper motion-matched 
point -6/ir = 0.05, the growth rate for the predicted motion matches that of the 
aerodynamics. The three motion-matched points are projected to a plot of p 
against V/bj„W 2 \/y. (The latter scale is at the upper right in fig. 3(b).) 

The question arises as to whether there is a significance to points on 
these curves that do not match the input 6. Such points can be interpreted to 
represent ordinate values of -6/tt + g, where g is the uniform-for-all-modes 
increment of g that would make those other points the ones of matched motion. 
For example, in figure 3(al on the curve segment for -6/ir = -0.05, the point 
with an ordinate of -0.03 represents (-6/tt) + g = -0.03, from which g = 0.02. 

Thus far only one value of k = 0.47 has been presented. Figure 4 pre- 
sents the results from solutions for a range of k-values from about 0.33 to 
about 0.55. The flutter-speed index V/bj^2Vv* plotted against density p. 
The three points from figure 2(b) obtained with simple harmonic aerodynamics 
(6 = 0) and the three points from figure 3(b) obtained with motion-matched aero- 
dynamics are indicated by the symbols. The point for 6=0 is common to the 
two sets. Solutions from a range of k give the isodecrement curves shown, 
solid-line curves for the three motion-matched values and dashed-line curves 
from a V-g type prediction that used simple harmonic aerodynamics. Curves 
for other decrement rates 6 are not plotted. 

Use of isodecrement results .- An example of the use of isodecrement curves 
is shown in figure 5. The isodecrement curves from figure 4 are included in 
figure 5. A possible wind-tunnel operating curve is shown in figure 13 of ref- 
erence 15. In this reference the vertical lines to the flutter boundary are for 
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fixed M and increasing dynamic pressure. An example of that general type of 
tunnel operating curve (fixed M and c and increasing P) is included in 
figure 5 rising from the lower left. The intersections of the operating curve 
with the isodecrement curves give the associated values of flutter-speed index 
and p. Projections can be made as shown in figure 5 to obtain the associated 
plot of decrement against density. On this plot the rate of change of the decre- 
ment with increasing p predicted from the matched aerodynamics (the solid-line 
curves) is readily seen to be sharply less (about one-third) than is predicted 
from V-g type curves with simple harmonic aerodynamics (the dashed-line curves) . 

The preceding results are for M = 1.3. Similar results have been calcu- 
lated over a moderate range of supersonic M and are shown in figure 6. Bound- 
aries of flutter-speed index are plotted against M. The contrast between the 
two types of prediction is virtually the same throughout the Mach number range. 


All-Movable Horizontal-Tail Flutter Model 

Flutter model HT-7 of figures 6 and 7 of reference 17 was also analyzed. 

It is a semispan model of an all-movable horizontal tail. In figure 7(a) the 
flutter-speed index is plotted against M, and the same type of boundaries as 
before are shown. The density P is 1.70 kg/m^, which is the density of the 
experiment at M = 1.64. For this case the contrast between the two types of 
prediction is even stronger; that is, the dashed-line curves from a V-g type 
prediction with simple harmonic aerodynamics (input 6=0) are very close to 
the flutter boundary, thereby indicating a steep gradient of motion decrement 
with changing M. (Although not presented here, the gradient-of-motion decre- 
ment with changing density for fixed M and c is also very steep from a V-g 
type prediction.) Observation of the solid-line curves from motion-matched 
aerodynamics shows a small rate of change of motion decrement with changing M. 
That rate of change is also small with changing P for fixed M and c. 

Figure 7 (b) gives the same results for model HT-7 in terms of stiffness- 
altitude parameter plotted against M as in reference 17. In figures 7(a) and 
7 (b) the upper Mach number limit of the calculated boundaries is the value of M 
for a sonic leading edge, which is the upper limit for the supersonic kernel- 
function ccmputer program used. 


C»NCLUDING REMARKS 

The equations of dynamic equilibrium are reviewed for a lifting surface in 
an air stream undergoing general exponential motion, that is, growing and decay- 
ing oscillations. With the appropriate motion-matched aerodynamics in these 
equations, analyses are made of preflutter (subcritical) motion and postflutter 
motion, as well as the usual flutter boundary. 

The three-dimensional supersonic velocity potential and kernel function 
were generalized to complex reduced frequencies appropriate for growing and 
decaying oscillatory motions. An existing computer program for the kernel- 
function method was modified for the eigensolution of the equilibrium equations. 
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Illustrative calculations were made for two flutter models for which these 
aerodynamic forces are applicable. In addition to the usual flutter boundary 
for zero motion decrement, example isodecrement curves were calculated for 
selected decaying and growing motion decrements. For these two models the rates 
of change of the motion decrement as functions of Mach number and air density 
were found to be much milder than the rates of change predicted from V-g type 
(velocity plotted against damping required) curves using simple harmonic (zero 
decrement) aerodynamics. 

This type of analysis has applications to the prediction of motion decre- 
ments during wind-tunnel and flight flutter testing and to the design of auto- 
matic controls for flutter suppression. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 31, 1978 
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APPENDIX 


THE RELATION OF SOLID-FRICTICMJ (MECHANICAL HYSTERESIS) DAMPING 

TO THE DAMPING DBCREMEMT OF EXPONENTIAL MOTION 

The differential equation of motion of a simple spring-mass-damper system 
subjected to a continuous simple harmonic forcing function is, in complex 
notation 


mx + cx + kx = 


(AT) 


where x is displacement, m, c, and k are constants representing mass, 
damper constant, and spring constant, respectively, F is the real amplitude 
of the forcing function, and a is the phase angle between F and x. The 
solution of equation (A1 ) is 

X = n - ei (wt+o-n) (A2) 

k 


where the magnification factor n is a function of m, c, and k, and n is 
the phase delay angle. 


For free oscillations (i.e., F = 0) , the logarithmic decrement 6 
ratio of amplitudes at the beginning and end of one ccanplete cycle is 


of the 


c/ccr c 

— 97f + n 


\jl - (c/Ccr)^ 

_\Ccr/ _ 


c 

2ir 

Ccr 


(A3) 


where c^j- is the critical value of c. 

For solid-friction structural damping (mechanical hysteresis) experimental 
results (see refs. 18 and 19) have shown that for forced simple harmonic motion 
the energy dissipated (converted to heat) by the structural damping is propor- 
tional to the motion amplitude squared but varies little with frequency. Conse- 
quently a good approximation is to make the replacement 

CO) = kg (A4) 

where g is the structural damping coefficient, and to treat k, g, and kg 
as constant. Thus for solid-friction damping 

mx + — i + kx = Fei(“t-*«) (A5) 

0) 
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for which 


X + Wn2(i 


+ ig)x = 


F 

_ ei (Wt-KX) 
k 


(A6) 


where 


0) 


n 


2 


k 


m 


The solution of equation (A6) has the same form as equation (A2) . 

For the case of free oscillations with solid-friction damping, the method 
described by Bishop in reference 20 is based on the assumption that equation (A4) 
is still an adequate approximation for exponentially damped oscillations, even 
though the motion is not simple harmonic. (In actual practice the common method 
of experimentally determining g is to measure the decay of a free oscillation.) 
Thus, 


6 = ■ = TTg + 0(g2) ~ TTg (A7) 

1 + \/l - g2 

This result is that of reference 21 and, through the order of g, is the same 
as that of references 22, 23, and others. 
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Figure 2,- Flutter and near-flutter eigensolutions obtained with simple harmonic 
aerodynamics (6 = 0) for clipped delta-wing model of reference 15. M = 1.3, 
k = 0.47, and gj^ = 0. Only critical (mode 2) root is shown. 
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-6/it plotted against p, and matched point for 6=0 and intersection 
points for -6/it = ±0.05 in plane of V/bjnO)2\/lJ plotted against p (at 
upper right) . 


Figure 2.- Concluded. 




(a) -6 /it plotted against 

Figure 3.- Near-flutter eigensolution curve segments obtained with motion- 
matched aerodynamics, input -6/Tr = ±0. 05, for clipped delta-wing model 
of reference 15, and 6=0 results from figure 2. M = 1.3 and 
k = 0.47. Only critical (mode 2) root is shown. 
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6/ir plotted against p, and matched points of V/bj/^2^ plotted 
against p (at upper right) . 


Figure 3.- Concluded. 
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Figure 4.- For clipped delta-wing model, flutter boundary, 6=0, and isodecre- 
ment curves, -6/tt = ±0.05, in plane of flutter-speed index plotted against 
density for M = 1.3 and range of k-values. Points for k = 0.47 from fig- 
ures 2(b) and 3(b) are indicated. 
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Figure 5.- Possible wind-tunnel operating curve through curves from figure 
and associated motion decrement plotted as function of density. 
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Figure 6.- For clipped delta-wing model, flutter and isodecrement curves in 
plane of flutter-speed index plotted against M for p = 1.225 kg/m^. 
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